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Abstract 

We demonstrate linear stability of the dilatonic Black Holes appearing in a 
string-inspired higher-derivative gravity theory with a Gauss-Bonnet curvature- 
squared term. The proof is accomplished by mapping the system to a one- 
dimensional Schrodinger problem which admits no bound states. This result is 
important in that it constitutes a linearly stable example of a black hole that 
bypasses the 'no-hair conjecture'. However, the dilaton hair is secondary in the 
sense that is is not accompanied by any new quantum number for the black hole 
solution. 
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1 Introduction 



In ref. [ l] we have presented analytic arguments in favour of, and demon- 
strated numerically, the existence of dilatonic black hole solutions with non- 
trivial scalar hair, in a string-inspired higher-derivative gravity theory with a 
Gauss Bonnet (GB) curvature-squared term. The numerical solutions clearly 
demonstrated the existence of a regular event horizon and asymptotic flat- 
ness of the four-dimensional spherically-symmetric space-time configurations 
considered in the analysis. There is a non-trivial dilaton (global) charge, 
which however is related to the ADM mass of the black hole, and hence 
the hair is of secondary type [0 : the gravitational field acts as a source for 
the scalar hair and one does not obtain a new independent set of quantum 
numbers characterizing the black holeQ. Subsequent to the work of Q , other 
researchers confirmed these results by discussing the internal structure of the 
solutions behind the horizon, and demonstrating numerically the existence 
of curvature singularities H . Also an extension of the analysis of ref. Jl|] to 
incorporate gauge fields became possible || [5J. 

The key feature for the existence of such hairy black holes is the bypass 
of the no-scalar-hair theorems || due to the fact that, as a result of the 
GB term, the scalar field stress tensor becomes negative near the horizon M, 
thereby violating one of the main assumptions in the proof of the no-hair the- 
orem |§. An additional element was the fact that the higher-derivative GB 
term provides a sort of 'repulsion' that balances the gravitational attraction 
of the standard Einstein terms, and a black hole is formed. In this respect 
the GB term plays a similar role to the non-Abelian gauge field kinetic terms 
in Einstein- Yang-Mills-Higgs theories which are also notable exceptions 
of the no-scalar- (Higgs)-hair theorem. 

An important question which arises concerns the stability of the dilatonic 
black holes. In the Yang-Mills case, the structures are similar to the sphaleron 
solutions of flat-space Yang-Mills theories, and thus unstable [Q . This may be 
easily understood by the fact that the black hole solutions owe their existence 
to a delicate balance between the gravitational attraction and the Yang-Mills 
repulsive forces. On the other hand, the dilatonic black hole solutions are 
entirely due to the existence of a single force, that of gravity. This already 

1 A similar situation occurs in the Higgs hair of the Einstein- Yang-Mills-Higgs sys- 
tems H D , where the Yang-Mills field act as a source for the non-trivial configurations of 
the Higgs field outside the horizon. 
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prompts one to think that such structures might be stable. 

It is the purpose of this article to argue that the dilatonic black holes 
are indeed stable under linear time-dependent perturbations of the classical 
solutions. To this end, we shall map the system of gravitational-dilaton equa- 
tions for spherically-symmetric solutions into a one- dimensional Schrodinger 
problem, where the instabilities are equivalent to bound states. We shall 
prove that our dilaton-graviton system admits no bound states. This result 
is important, since it constitutes an example of a hairy black hole structure 
that appears to be, at least linearly, stable. Its importance is also related 
to the fact that such higher-curvature gravity theories are effective theories 
obtained from superstrings, which may imply that there is plenty of room 
in the gravitational sector of string theory to allow for physically sensible 
situations that are not covered by the no-hair theorem as stated ||. Unfor- 
tunately, at present non-linear stability of the dilaton-graviton-GB system 
cannot be checked analytically, and is left for future investigations. 



2 Relevant Formalism 



We start by considering the action of the Einstein-Dilaton-Gauss-Bonnet 
(EDGB) theory: 

'R a'e* 
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where 
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GB 



The spherically-symmetric ansatz for the metric takes the form 

ds 2 = e r ^dt 2 - e A ^dr 2 - r 2 (d6 2 + sin 2 6 dip 2 ) 
The equations of motion derived from ([]]) are: 
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For later use we note that for the derivatives of the dilaton field at the horizon 
in the static case n one has the following behaviour: 
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We now notice that M only one of the two branches of solutions in (|9|), the 
one with the + sign, leads to asymptotic flatness of the fields, and this is the 
branch we shall consider here. The equation for <p" near is: 

<j>" = --^—L — — - r' + o(i) (n) 

2 -6^e 2 ^ + ^e^V 3 + 2r 4 

which is finite 0(1), as a result of (|9|). The asymptotic form of the dilaton 
field and the metric components near the even horizon r ~ are: 

e" A M = A 1 (r-r,) + A 2 (r-r,) 2 + ... 
e r(r) = 71 (r -r h ) +7 2 (r - r ft ) 2 + ... 

0(r) = /l + 0' /i (r-r,,) + , / :(r-r ?i ) 2 + ... (12) 



where: 



X 1 = 2/(a'e^4>' h /g 2 + 2r h ) , (13) 



and 71 is an arbitrary finite positive integration constant, which cannot be 
fixed by the equations of motion, since the latter involve only T'(r) and 
not T(r). This constant is fixed by the asymptotic limit of the solutions at 
infinity. At infinity, one uses the following asymptotic behaviour: 

iy \ /yd 

m - (u) 

which guarantees asymptotic flatness of the space time. Above, M denotes 
the ADM mass of the black hole, and D the dilaton charge. The numerical 
analysis of has shown that M and D are not independent quantities, 
thereby leading to the secondary nature of the dilaton hair (l|. [2|] . 

The black hole solutions of are characterized uniquely by two param- 
eters (4>h, r h)- Note, however, that the equations of motion remain invariant 
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under a shift — > + 0o as long as it is accompanied by a radial rescaling 
r — > re^°/ 2 . Due to the above invariance it is sufficient to vary only one of 
and 4>h- In the present analysis we choose to keep r h fixed {r h = 1), and to 
vary A typical family of solutions has dilaton configurations (outside the 
horizon) of the form depicted in Figure 1. The solutions are characterized 
by negative <ph, and monotonic, non-intersecting behaviour from r/j until in- 
finity. These are essential features of the solutions, that we shall make use 
of in our linear stability analysis. 

3 Linear Stability Analysis 

We now consider perturbing the equations ([|)-(|D by time-dependent linear 
perturbations of the form: 

T(r,t) = r(r) + ST(r,t) = T(r) + 5T(r)e iat 

A(r, t) = A(r) + 8A(r, t) = A(r) + 5A(r)e iat 

(j)(r,t) = <t>(r) + 6<t>(r,t) = (j)(r) + 5<j)(r)e iat (15) 

where the variations ST, 5 A and S(p are assumed small (bounded), and the 
quantities without a 5 prefactor denote classical time-independent solutions 
of the equations (§)-(§)• The above harmonic time dependence is sufficient 
for a linear stability analysis, since by assumption the linear variations are 
characterized by a well-defined Fourier expansion in time t f|, |9|. The linear 
stability analysis proceeds by mapping the algebraic system of variations of 
the equations of motion under consideration to a stationary one-dimensional 
Schrodinger problem, in an appropriate potential well, in which the 'squared 
frequencies' cr 2 will constitute the energy eigenvalues. In the present problem, 
the 'wavef unction' turns out to be the dilaton linear variation 8(j){r)e lat . In- 
stabilities, then, correspond to negative energy eigenstates ('bound states'), 
i.e. imaginary frequencies a. As we shall show in this article, for the sys- 
tem of variations corresponding to (0)-(@),(0) the corresponding Schrodinger 
problem admits no bound states, thereby proving the linear stability of the 
EDGB hairy black holes. 

As we shall discuss below, some technical complications arise, as usual 0], 
in the above process due to the fact that the 'naive' stationary Schrodinger 
equation with respect to the original coordinate r is not well defined at some 
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points of the domain of r e [r^, oo). This necessitates a change of coordinates 
in such a way that the resulting Schrodinger problem is well defined. A 
convenient choice is provided by the so-called 'tortoise' coordinate r* j|, [J, 
which is defined in such a way so that the domain [r^, oo) is extended over 
the entire real axis r — > r* G (—00,00). In our specific problem, we shall 
define the 'tortoise' co-ordinate r* as H: 



dr* 
dr 



-Cr-A)/2 



(16) 



+ V(r*)u(r*) = — (T 2 -u(r*) 



(17) 



As we shall show, then, the associated stationary Schrodinger equation, per- 
taining to the dilaton variation in fll5|) , will be of the form: 

_d_ 

dr* 

where u(r*) is related to the dilaton variation 5<p(r) in (15), and the potential 
V is well defined over the domain of validity of r*. 

Let us now proceed with our analysis. The perturbed equations 
under the variations (|15[) read: 
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We can integrate eq. (|2lD to obtain 
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where /x(r) is an arbitrary function of r which can be set equal to zero if we 
require that SA = when 5(p = 0. An independent check of this assertion 
can be obtained as follows: first we differentiate (B^) with respect to r and 
then use eg. (|T9|) , as well as the time- independent equations of motion. Thus, 
we obtain the following differential equation for /x(r) 

/i'(r) + M(r) (y - y + I) = (24) 

This can be integrated to give /i(r) ~ e( A ~ r " 2 /r. When r — > r^, e^ A ~ r ^ 2 — ► 
oo and fi — > oo, which is incompatible with the assumption of a small <5A(r) 
required for the linear stability analysis. Rejecting this solution, we are 
left with the trivial one, \x = 0. For calculational convenience we also set 
a'/g 2 = 1, from now on, which is achieved by an appropriate rescaling of the 
dilaton field. 

Rearranging the above equations, one obtains, after a tedious but straight- 
forward procedure, an equation for 5(f) which has the following structure: 

A5(f)" + 2B5(j)' + C5<f) + a 2 E5(j) = (25) 

where A, B, C, and E are rather complicated functions of 0, 0', 0", A, A', 
T, r" and T". In the limit r — > these coefficients take the form 
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where we have used the asymptotic behaviour (|T2| ) near the event horizon. 
On the other hand, when r — > oo we obtain 
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where we have used the asymptotic behaviour fll4|) near infinity. 

As we can see from (p6|)-(p9|) the coefficients of the Schrodinger equation 
( |25|) are not finite at the boundary r = r^, where the variation <50 is bounded. 
As mentioned previously, to arrive at a well-defined Schrodinger problem, one 
can use the 'tortoise' coordinate (|T6|). Then, the perturbed equation for the 
dilaton field takes the form 
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where Ai is given by ([13]). As a result, all the coefficients in eq. ([34]) are now 
well-behaved near the horizon r^. In order to eliminate the term proportional 
to S(j)', we first divide eq. (B3) by A and then we use the function 




Then, the equation for 



= exp 

takes the form 
.£ B 2 



(40) 
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and we have set u = F 5<j). 

It is straightforward to see that 
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independently of r^, 0^. In addition 
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Moreover, as shown in Fig. 2, the function B/A is well-behaved over the 
entire domain outside the horizon, implying the integrability of the function 
F . Also, the quantity 8/ A is finite and always positive outside the horizon of 
the numerical black- hole solutions of ]IJ (see Fig. 3). It is also immediately 
seen that the eigenfunction Uq, corresponding to the eigenvalue a = 0, can 
be constructed out of the difference of any two curves in Fig. 1. From the 
monotonic and non-intersecting nature of the various members of the family 
of the numerical solutions of Fig. 1, then, one can conclude that Uq has no 
nodes in the domain r* G (— oo, oo), and that pIuq/uq = e r_A [|(r' — A')u' + 
u'q]/uq is finite. This, together with the finite and smooth form of B/A 
(Fig. 2), implies, on account of (41), the finiteness of the coefficient C/A 
outside the horizon, without the need for an explicit numerical computation. 



Thus, equation ([41]) assumes the form of an ordinary Schrodinger with regular 
coefficients over the entire domain of r*. 
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From (|4"0D, ( f43|) one obtains on the boundaries: 

P*Mo|r*=±oo = -^U \ r * =±oo + FpJ(f)\ r * =±00 (45) 

The 5<f> remains bounded at r* = ±00. From the asymptotic behaviour 
( |j~4"D it becomes clear that, at the r = 00 boundary, B/A — > B/A ~ 1/r, 

independently of <f>h- Thus, F(r* — > 00) = e^-°° A r ~ r for r — > 00. Hence, 
the first term in ( f45|) becomes 50oo + O(-), whilst the second term vanishes 
as 1/r, for r — > 00 (see fll4]) ). Hence, at the r = 00 boundary the boundary 
values of p*u are proportional to those of u : 

1 - 2£2.i ^ 

P*Uq r*=oo r*=oo 1 4tOJ 

r 




On the other hand, from fl40D and (J4j) it becomes clear that at the boundary 
r* = —00 (horizon), F{r* = —00) = 1. Moreover, p*Wo|r*=-oo 

p*«o|r*=-oo = (r - r h )(0 (2) ' - (1) ')| r=rh = ( (r - r h )^^u 

\ v<Ph 

since, in our construction, each family of solutions is uniquely character- 
ized |J by the value <f>h for fixed (we remind the reader that here we chose 
to keep Th = 1 (fixed) and vary 0^). From (|9|) one easily observes that for 
linear variations, 5cf> = 0^ 2 - ) — (f>^\ the difference 0^ — <p£ ls finite. Hence, 

p*uo|r*=-oo -> (48) 

We now discuss the a 2 < unstable modes iv It can be easily seen from 
(|30|) - (p3|) that the asymptotic form of equation (|4]) at the r* = 00 boundary 
reads: 

p\u a = -a 2 u a (49) 
Hence, the bound-state solution u a behaves as Q 

u a {r* = 00) = e- |<T|r * -> (50) 

2 Here, and in the following, wc insist on bounded, or- at most- linearly divergent u a , 
at the boundaries r* — ±00. This is due to the fact that, since u = FS<fi, and F is 
independent of a and at most linearly divergent at r* = 00, then it is only for such 
a behaviour that the variation 8<f> remains bounded, as required by the linear stability 
analysis. An exponentially divergent u„ ~ e r at the boundaries, would imply 8<f) ~ e r / F, 
and, hence, is not acceptable. 
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On the other hand, from (|26|)-(p9|) it is obvious that, on the horizon, equation 
(f4TD assumes the form (for the case = 1): 

plu a + k 2 u a = (51) 

k 2 = ; 716 + a 2 = kl + a 2 (52) 

(1 + Vl-Qe 2 ^)Vl-Qe 2 ^ 

From (|5^) one can see two possibilities near the horizon of the black hole: 

• (i) The 'total energy' is such that > a 2 > , 

Taking into account the asymptotic form at r* = oo ( pE§D one also 
observes that in this case the spectrum of the respective Schrodinger 
equation is continuous and non degenerate. The general solution of the 
perturbation u a near the horizon is, therefore, oscillatory (unbound 
state): 

Such a continuum of states cannot exist in our case by continuity. In- 
deed, due to the non-degenerate nature of the eigenvalue problem, the 
limiting case a —>■ should yield the solution u . However, in the limit 
a 2 — > 0, and in terms of r, one obtains 

/ e 4>n \ 

U ° ~ C ° S ^ (l-6e 2 ^) 1/4 ~ ^ + if °) 

where ip is a constant phase shift, and we have taken u to be the real 
part of eq. (|53D. Then, 



p*u k sin ln(r - r h ) + tp Q (55) 




The above result is not in agreement with eq. (f±8D, thereby contra- 
dicting the non-degenerate nature of these solutions, which, in turn, 
implies the absense of such solutions in the problem (|4l|). 

(ii) This leaves one with the second possibility of a discrete spectrum 
of bound states, which would occur for: 

a 2 < 2lie24 ,h ^ = < (56) 

(1 + VI - 6e 2 ^)Vl - Qe 2 ^ 
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As we shall show below this is also not realized due to the special form 
of Uq. 



To this end, one first observes that such bound states would vanish expo- 
nentially at the r* = — oo boundary: 



-oo) ~ e |fc|r * -> (57) 



Thus, on account of (|||f|8[|50],|57|) the Wronskian of any two solutions u±,u 2 



of the equation fl4l|) with a 2 < vanishes at the boundaries: 

W = {UiP*U 2 - M 2 P*Ml)|r*=±oo = (58) 

To count the unstable gravitational modes of the original problem, one 
needs to count the nodes of the wave function u of the one-dimensional 
Schrodinger problem ([O]). Fortunately, this can be done without detailed 
knowledge of the solutions. As we shall discuss below, all one needs to 
observe is the monotonic and non-intersecting nature of the dilaton curves 
in Fig. 1. To this end, one first observes that a standard 'node rule' for the 
discrete spectrum of the equation (|4l|) applies, which is a direct consequence 
of Fubini's theorem of ordinary differential equations [TO]. This theorem can 



be stated as follows: consider two differential equations: 

u" + 2 Pl u' + q lU = (59) 
u" + 2p 2 u' + q 2 u = (60) 

If, 

P 2 +P 2 2 -Q2< A +Px~qi (61) 

throughout the interval [a, b], then, between any two consecutive zeroes of a 
solution of fl59D, in the interval [a, b], there is at least one zero of a solution 

In our case, we can apply this theorem for two-different eigenfunctions, 
Ui, u 2 , corresponding to eigenvalues a\ and of of ([4lD. The interval [a, b] is 
the entire domain of validity of the solutions of (|4T|), (— oo, oo), including the 
boundaries at infinity. In this case, 

P, = 0, qi =^ + a^-^-p*(^), i = l,2 (62) 
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Then, the (sufficient) condition for a 'node rule' (|6T| ) reads simply: 

£ o £ 



-07, 



a^-T' (63) 

The positivity of £/A (Fig. 3) implies that the condition fl63|) becomes 
simply: 

of > a? (64) 

This special version of the theorem is known as Sturm's theorem ||1 1|| . As 
a corollary of Fubini/Sturm's theorem one obtains the standard 'node rule' 
for the number of zeroes of the eigenfunctions in the discrete spectrum of 
bound states, according to which if the eigenfunctions are ranked in order of 
increasing energy, then, the n — th eigenfunction has n — 1 nodes (excluding 
the boundary zeroes) [[11] ('node rule'). 

Consider, now, the case where 02 corresponds to the zero eigenvalue of 
(f4i"D, <7o = 0. As can be seen from the numerical solution of Fig. 1, the 
monotonic and non-intersecting character of the dilaton curves in the entire 
domain outside the horizon implies that the solution w 2 = m , which, as we 
have mentioned earlier, can be constructed out of the difference of two such 
solutions, has no nodes in the domain r* G (—00, 00). Since any solution u n 
from the discrete set of negative eigenvalues a 2 < (unstable modes) has 
at least two nodes at the boundaries, according to Fubini/Sturm's theorem, 
«o should have at least one in the domain r* G (—00,00). This contradicts 
the fact that u is nodeless. Thus, the only consistent situation is the one 
without such negative- energy modes. This, in turn, implies linear stability 
for the dilaton- GB black holes of ref. [p]]. The reader might worry about 
the divergent boundary conditions of «o at r* = 00, which makes it not 
an ordinary eigenfunction of a Schrodinger problem. In the appendix we 
argue that this is not an obstacle. In fact, we present an explicit proof of 
the absence of bound states in our case, following the same spirit used in 
the proof of Fubini/Sturm's theorem. The crucial element, which allows 
the standard proof to go through, is the special boundary condition for the 
Wronskian (pB[), which is valid for the entire spectrum of eigenfunctions of 
( f4T|) with o 2 < 0, including the non-standard one uq. 

The above considerations can be extended straightforwardly to the case 
where r^ — > 0. All the coefficients of ( f4T|) are still well-defined in this case, 
which implies that the stability in principle does not change. However, ac- 
cording to the analysis of ref. M, the case r/j — > corresponds to a singular 



14 



curvature scalar 

R ~ \ r h -> (65) 

which implies the absence of 'particle-like' solutions in the EDGB system |E], 
From the condition for the existence of black hole solutions (TO), one 



observes that the only consistent value of ^ for r fc — > is [J] 4>h — > — oo. 
The limit is taken in such a way that fta/ey g 2 r% = 1. This implies that in 
such the GB term in the action ([I]) becomes irrelevant, and one is 

left with the standard Einstein term, which admits only Schwarzschild black 
holes, known to be stable. This stability is confirmed by the above smooth 
limit of the coefficients in (|4lD as r/, 0. 



4 Remarks and Outlook 

Above we have demonstrated linear stability of the dilatonic black hole solu- 
tions in the EDGB system, found in ref. 0. This result is important, since 
it constitutes an example of a stable, albeit secondary, hair that bypasses the 
no-hair conjecture 0, |8[. Non-linear stability of the EDGB system, however, 
although expected, still remains an open issue. 

Before closing we would like to compare our semi-analytic results on linear 
stability with some remarks in favour of stability by virtue of a catastrophe 
theory approach made in ref. ||. As usual |IBfl , catastrophe theory can 
only indicate relative changes of stability, and hence cannot constitute a 
'proof of stability. In ref. || a numerical solution was found for the (+)- 
branch of solutions (0) which indicated the existence of a 'turning point' (TP) 
in the r^-M (or equivalently 0^-M) graph. The TP occurs at the 'critical 
point' for the existence of black hole solutions, which is the point where the 
black hole acquires a minimal mass, below which no solution is found. In 
the numerical solution of ref. |J a continuation beyond this critical point 
emerged, which ends at a point ('singular point') where a singularity appears 
in the square of the Riemann tensor, as well as in §" h . 

The part of the solution from the critical point to the singular point was 
argued in ref. || to be relatively unstable, as compared to the regular branch 
discussed here. Such a change in stability manifests itself as a cusp in an 
appropriate catastrophe theory diagram [|| [13|]. In ref. || such a diagram 
has been chosen to be the diagram of the thermodynamic entropy |14| versus 
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the mass of the black hole. 

In our numerical solutions |l[], used here, we found no evidence for such a 
continuation of the solution after the critical point, of minimum mass, which 
in our analysis is also the end point, where the two branches of (|9|) meet (see 
Figure 3). Our black hole solutions are uniquely specified by the pair (r h , 4> h ), 
which was essential in our linear stability analysis above. Moreover, for finite 
rh, the quantity §" h is finite In this respect we are in agreement with the 
results of ref. 0, where a branching of solutions was found only inside 
the event horizon. These authors have also given a graph of the entropy 
versus the mass of the black hole solution of ]I[ outside the horizon, and 
found, as expected, a smooth curve, with no cusps. In this article we have 
proven analytically the stability under linear perturbations of this (unique) 
branch of the black hole solutions, which in the — M graphs of Figs. 
4,5 terminates at the minimum-mass critical point. Of course, this result 
is in agreement with the relative stability of this branch suggested by the 
catastrophe theory analysis of ref. |§, but, as we said, the form of their 
numerical solution appears to be different, as it continues beyond the critical 
point. Nevertheless, we consider this as an indication for the stability of the 
solutions beyond the linear approximation. This, however, still remains an 
open issue. 

As a final remark we would like to mention that the effects of gauge 
fields on the dilatonic black hole GB solutions have been considered in ref. 
[§, 0]. From the point of view of stability, one expects that, in the case of 
'coloured' black holes, involving non-Abelian gauge fields, instabilities occur 
in both the gauge and gravitational sectors of the solutions. Instabilities 
in the gauge sector are of sphaleron type [|J. Those in the gravitational 
sector can be studied in a similar way as for coloured black holes in Einstein- 
Yang-Mills theories H. One can go beyond linear stability analysis in such 



systems, by invoking catastrophe theory [13| , which is capable of giving 
the relative stability of various branches of solutions for the coloured EDGB 
black holes || [7| . However, analytic methods can still be combined with the 
catastrophe theory approach || in order to count the unstable modes in both 
sectors, gauge and gravitational, by invoking appropriate maps of the system 
of perturbations into one-dimensional stationary Schrodinger problems @. [|. 
We hope to return to a detailed analytic study of these issues in a future 
publication. 
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Appendix: Absence of bound states 

In this Appendix we prove the absence of bound states in the problem fl4lD by 
following a Wronskian treatment for the entire set of eigenmodes with a 2 < 0. 
This justifies the validity of Fubini/Sturm's theorem in our case. We stress 
that the crucial point in the proof is the special boundary conditions of the 
Wronskian (|58|) . These allow a standard Wronskian treatment to go through 
for the u solution of fl¥Ip, despite its (linear) divergence at the r* = oo 
boundary, which makes it not an ordinary eigenfunction of a Schrodinger 
problem. 



To this end, we consider first the two solutions of fl4ip , tt , an d Ub - the 
ground state, at the bottom of the discrete spectrum. Both of these have 
no nodes in the interior domain of r*, excluding the boundaries (the node 
structure of Ub follows from the 'node rule'). We, then, employ properties 
of the Wronskian of the solutions as follows: first we multiply each equation 
with the other eigenfunction. Next, we subtract the resulting system of 
equations, and then, integrate it over the entire domain of r* e (—00,00). 



In this way one obtains, in a standard fashion |l I 



f°° £ 

AW\ r * =±OQ = {al - ct'q) I dr*—u b u (66) 



00 



where the left-hand-side denotes the change in the Wronskian between the 
two boundaries. From (|58|) this vanishes. Moreover, as we have mentioned 
previously, 8/ A is positive definite for the entire domain of r* G (—00,00) 
(Figure 3). Since Ub, u have no nodes in the domain (—00,00), excluding 
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the boundaries, one obtains from fl66|) that the only consistent case is the 
degenerate one of = Oq. But Oq = 0, whilst of < by assumption; this 
implies a contradiction, excluding o& from the spectrum. 

One repeats the construction, using u and any of the higher eigenfunc- 
tions of the discrete spectrum, u n , corresponding to a\ < 0. The change in 
the Wronskian between — oo and the first encountered zero of u n , at r* = zq, 
is then given by: 

AVT|^ = _ 00 = -u p*u n \ zo = {a 2 n - al) / dr*-u n u (67) 

Without loss of generality, one may assume that u n > in the interval 
(— 00,20) ■ Then p*u n (z ) < 0. It is immediate to see that there is a contra- 
diction in (|67D . The middle part has the sign of u , whilst the right-hand-side 
has the opposite sign of u . The case of a zero of u n , at a point z , such that 
p*u(z ) = is dealt with similarly. In that case the contradiction lies in the 
fact that the left-hand side vanishes, whilst the right-hand side is a negative 
number (for u n > in the interval). These results exclude the possibility of 
bound-state eigenfunctions with zeroes in (—00, 00). 

The above analysis, therefore, implies the absence of negative energy 
modes (bound states) in the problem fljip, which, in turn, leads to linear 
stability for the Dilaton-Gauss-Bonnet black holes of ref. |JJ. 
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Figure 1: Dilaton configurations for a family of black hole solutions, corre- 
sponding to fixed Th — 1, for various values of 4>h- Notice the monotonic 
behaviour of the solutions, and the fact that the curves do not intersect. 
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Figure 2: The graph depicts the coefficients B/A and B/A, for a typical 
member of the family of the black hole solutions of Fig.l corresponding to 
<f> h — — 1, r h — 1. It is clear that the coefficient B/A, incorporating the 
tortoise coordinate, is finite in the entire domain outside the horizon, thereby 
implying that the quantity F is well-defined and integrable. 



22 




1 10 10 2 10 3 

r 

Figure 3: This diagram depicts the coefficients E/A and E/A for a typical 
member of the family of the black hole solutions depicted in Fig. 1 ((p^ = — 1, 
r h = 1); the coefficient E/A diverges at the horizon as l/(r — r h ) 2 . On the 
other hand, S/A = e r ~ A E/A, appearing in (40), is finite at the horizon. The 
positive-definiteness of both coefficients is clear. 
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Figure 4: The graph depicts r h versus the ADM Mass 2M of the black hole, 
for a fixed value of 0^ = —1. The emergence of the asymptotic critical point, 
r\ ~ 6a 2 e 2< t >h /g 4 , below which there are no solutions, is apparent. At this 
point the mass becomes minimal. 
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Figure 5: The magnification around the critical point, depicted in the above 
figure, shows clearly the abrupt (almost vertical) slope with which this point 
is approached. No turning point is found. 



25 



